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Abstract 

The phase diagram and low temperature properties of the doped manganites 
Ai-ajBjjMnOs are discussed for the concentrations x < 0.4. The transition 
from insulating antiferromagnetic to metallic ferromagnetic state at x cr ~ 0.16 
is treated by means of percolation theory. The unifying description of insu- 
lating and metallic states is presented. The undoped manganite is a band 
insulator consisting of ferromagnetic layers, which are coupled antiferromag- 
netically along the c direction with a low Neel temperature. The metallic 
phase can be described by the two-band Fermi liquid picture. The behavior 
of conductivity, spin wave excitations, etc. is analyzed and the comparison 

with experimental data is carried out. 
PACS numbers: 72.15. Gd, 75.30.Ds 
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I. INTRODUCTION 



Interest in utilizing the "colossal magnetoresistance" (CMR) effect has resulted in enor- 
mous literature, which spans properties of manganites at numerous compositions, ceramic 
and crystals, crystalline films on different substrates, stoichiometric and not (for extensive 
review see Nevertheless, until recently, it was somewhat difficult to establish sys- 

tematic trends which govern the physical properties of these materials (especially their low 
temperature behavior), to access peculiarities pertinent to the ground state and its depen- 
dence on composition. 

As for the overall behavior, especially for the high temperature properties, which are 
responsible for the CMR phenomena itself, the consensus is that the Zener's double exchange 
(DE) mechanism together with the Jahn- Teller ( JT) effect for Mn 3+ as suggested in ||[7| 
produce the basis for their understanding. High temperature lattice deformations localize 
carriers by a formation of the JT or local polarons |6|||. With temperature decrease, DE 
mechanism delocalizes carriers and result in simultaneous onset of metallic conductivity and 
half-metallic ferromagnetic state. The physical picture, is by all means correct, but the 
method, introduced in ||, is not sufficient for understanding delicate qualitative features, 
especially at the low temperatures. 

As an important issue, we have to mention, is the evolution of ground state for the 
doped manganites from the insulating antiferromagnetic state of LaMnOs to the pronounced 
metallic ferromagnetic one in Lai-^Sr^MnOs at x ~ 0.3. 

It has been demonstrated first in our short paper || that both the half-metallic state at 
finite concentrations and the insulating state of the parent LaMnOs can be brought together 
by making use of a simple band model. Parent LaMnOs is considered in that scheme as 
a band insulator. As for the doping process, it has been shown in ||, that it leads to the 
percolative scenario, producing the well defined threshold concentration x cr ~ 0.16 — 0.17 
for the onset of low-temperature metallic conductivity. 

In the presentation below we describe in more details and further elaborate the theory, 
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proposed by Gor'kov and Kresin, keeping in mind both some future applications and anal- 
ysis of available experimental data at x < 0.4 is also provided. It is shown that the two 
bands theory agrees reasonably well with the results of recent experiments, demonstrating 
applicability of the Fermi liquid approach to the metallic manganites. The scheme serves 
as a theory unifying rather rich low temperature properties of manganites and is capable 
of some predictions regarding the evolution of physical parameters, the shape of the Fermi 
surface etc. with concentration. The experimental data show also the consistency with our 
results in the metallic concentration range and trace remaining percolation features. 

The paper is written as follows. In Sec. II the main features are identified and accounted 
in Hamiltonian. Sec. Ill contains the energy spectrum and its dependence on a spin ordering 
or JT - deformations. Properties of a stoichiometric compound LaMnO"3, treated as a 
band insulator, are considered in Sec. IV. Growth of the percolative regime for the phase 
separation and the critical composition for onset of metallic conduction are discussed in Sec. 
V. Sec. VI contains a summary of major metallic properties calculated in the framework 
of the two band model, while in Sec. VII a comparison with available experimental data is 
discussed. Sec. VIII contains brief summary and conclusions. 

II. MODEL AND HAMILTONIAN 

We just start with, a somewhat oversimplified (the so-called "pseudocubic" ) crystalline 
structure for manganites. As it is well-known in that approximation the unit cell, 

say, for LaMnC>3, may be taken as a cube, with the lattice constant ao — 3.9A Rare- 
earth or alkaline ions are placed at the center , while the manganese ions occupy the corner 
sites. Mn 3+ - sites are caged into the oxygen octahedra, which share the O 2- ions along the 
Mn-O-Mn bond (another view would be that each La 3+ - ion is cooped up in the midst of 
twelve O 2- ions). The ideal structure is then modified for real materials, AMJ1O3, due to a 
mismatch in the ionic radii. The latter is commonly characterized by the tolerance factor, t 
(see e.g. in 0): 
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t — — (1) 

V2 i?Afn + Ro 

The effect of t 7^ 1 is that the oxygen octahedra become periodically tilted, and the unit 
cell may then be comprised of a few "pseudocubic" cells. It is shown below that deviations 
from the "pseudocubic" structure, i.e. deviations in the angle, a, of the Mn-O-Mn bond 
from 180° are not of much importance for the "average" electronic structure. However, we 
will see later that local fluctuations in the tolerance factor ([!]) may play rather significant 
role for the conducting properties of the "doped" manganites Ai-^B^MnOs. 

Let us outline our base model which is then applied for the interpretation of manganites 
properties. Since the manganese d- shell has reasonably small orbitals, the model exploits 
the tight binding approximation. 

As it is well known, the Mn 3+ ion has four d- electrons. The d- shell in the cubic 
environment is split into the triple (£2$) and double (e 2g ) degenerate terms. The t 2g - level is 
fully occupied by three electrons considered as a local spin S = 3/2, in accordance with the 
Hund's rule. The e 2g - term in manganites may be empty (Mn 4+ ) or single occupied (Mn 3+ ). 
If there is one electron on the e 2g - level, the direction of its spin is governed by the same 
intra-atomic Hund's interaction: 

H H = -J h <t- (2) 

(<r- the Pauli matrix for the e 2g - electron). The Hund's coupling Jh > is rather large 
(JhS ~ 1 -T- l.beV). Therefore spins of e 2g and t 2g electrons are ferromagnetically aligned 
locally. 

The two major assumptions are: 1. - the splitting between the t 2g — and e 2g — levels 
is large enough, so that only e 2g electrons may participate in conductivity; 2. - the e 2g — 
levels are well above the top of the filled up oxygen band (on a scale of a few eV). The 
band-structure calculations do not always support the last assertion (see e.g. |ID|). On the 



other hand, extensive optical studies agree well with the second assumption, see e.g. [[LI]]. 

Hoping process of one electron from a site % to its nearest neighbor % + 5 in the tight- 
binding approximation has to be modified to account for the double degeneracy of the e 2g - 



level on each site. In the hopping Hamiltonian, hence: 



Ht — ti t i + s, 

i,S 



(3) 



i^i+s becomes a two-by-two matrix. Its explicit form depends on choice of the basis for the 
e2g - representation. 

The C2 g - degeneracy of the Mn 3+ site would locally result in the Jahn - Teller (JT) 
instability with respect to spontaneous lattice distortions. Two active modes for deformation 
of the surrounding oxygen octahedron are usually denoted in literature as Q2 and Q3 (see 



12 and review 111311), and form the two-dimensional basis of the octahedron deformations. 



In the invariant form the JT part of the electron-lattice interaction may be written as: 



Hjt = 9 fi • Qi 



(4) 



where f is called the "pseudospin" matrix (see fll3l ). 

In what follows we use the normalized complex basis for the two-dimensional 

representation e 2g : 



ipi oc z 2 + ex 2 + e 2 y 2 ; i/j 2 = 



(5) 



where e = exp(27rz/3). Another basis, which is often used in literature, for example P~5| , |r5 
is given by the real functions: 



if 1 oc d z 2 = 3z 2 — (x 2 + y 2 ); tp 2 oc d x 2_ y 2 = x 2 — y 2 



(5') 



Connection between these two is: 



V>i = {fx + if 2) tp 2 = 4>l- 

The given basis (^j) makes it possible to present (^) in the convenient form: 

( n \ 



(5") 



exp(il 
exp(-iO) 



(6) 



/ 



where in the standard notations | E||13Ul5 | 



Q2 = Qo sin j Q3 = Qo cos 6 (7) 

with Qo being the magnitude of the JT - distortions. The angle, 0, specifies the shape of 
the distorted octahedron. Thus, the angles 9 = 0, ±27r/3 correspond to elongation of the 
octahedron along the z, x and y - axes, respectively. 

The JT - term (|]) is linear in Q, while the elastic energy is quadratic in Q. Therefore if 
one electron is placed on the JT - level, the total site energy always decreases with non-zero 
lattice deformations. Strictly speaking, deformations, Qi, on the two adjacent manganese 
sites are not independent because the two sites share one oxygen along the Mn-O-Mn bond; 
in what follows only cooperative JT distortions will be considered. 

According to review [JT2|| (see also we are dealing with two normal modes Q 2 and 



Q3, where: 

Q2 = -75O1 - x 4 + y$ - y 2 ), 

(8) 



^ 3 = 7E( 2z2 ~ 2ze ~ ^1- x 4 -y 2 + y 5 )- 
Onset of the cooperative JT distortions means the onset of a non-cubic structural order 
( "ferrodistorsive" deformations, if the wave vector of the structure, q = 0, "antiferrodistor- 
sive" ones, if q 7^ |12|,|16j). To account for fluctuations or changes in the phonon spectrum, 
the displacements in (|5p must be expressed in terms of the normal modes for the crystal 
vibrations. 

With the above remark regarding the JT-distortion, one may now bring all the contri- 
butions (HHH) together: 



(9) 



(the last term added into @ is responsible for the "elastic" energy of the JT mode). The 
Hamiltonian of exactly the form (Q) has been considered by Millis et. al. As it 

was mentioned above, the important physical observation made in [||,[7| was that with the 
temperature increase disorder produced by thermally excited JT degrees of freedom tends 
to localize the charge carriers. Though the idea by all means is qualitatively correct, the 
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approach, used in |J, is a selfconsistent interpolative scheme which is not appropriate for 
understanding low temperature properties of manganites, i. e. the symmetry of the ground 
state and its variation with composition. The latter issues are addressed below mainly in 
terms of the band approach. 



III. BAND SPECTRUM 



Let us consider first the band spectrum of the Hamiltonian (the Hund's term and 
JT term being temporary omitted). The matrix i in (|3|) on the basis (|5|) has the form: 



^11 ^12 



J 21 ^22 



(10) 



where 

S n = S 2 2 = S = (A + B)[cos(k x a) + cos(k y a) + cos(fc^a)], 
E12 = = (A — B)[cos(k z a) + ecos(k x a) + e 2 cos(k y a)}. 

A and B in ([11]) are two overlap integrals: 



(11) 



A oc tpxfa x, y)ifi(z + a; x, y), 
B oc ip 2 (z; x, y)ip 2 (z + a; x, y), 



(12) 



where the bar (...) means the matrix elements for the interaction potential on the two 
Wannier functions, <fi,2, of the neighboring atoms (in the tight binding approximation A = 
— \A\ < 0). Simple geometric considerations for the d-shell show that \B\ < \A\(\B\ ~ ^|^4|, 
according to ||T7l , |T8l1 ). The cubic spectrum consists of the two branches: 



£i,2(p) = + |S|)(cb. + c y + c z ) ± {\A\ - \B\)^c 2 x + c 2 y + c 2 z - c x c y - c y c z - c z c x (13) 

(we introduced the notations q = cos(fcja), % = x, y, z). 

If a homogeneous "ferroelastic" JT-distortion is imposed on the lattice, the spectrum 
may be obtained by adding —SRi eX p term to S i2 in (|TUp. For some future applications 
we provide here the spectrum for the limiting case \gQo\ ^> \A\: 



£!, 2 (p) = ± | ^ | +(A+ B)(c x + c y + c z ) ± 

(A - B) j(c z - \{c x + c,)) cos (B) + ^§(c y - c x ) sin (61) | . (14) 

This expression is interesting: at the choice sin (6) = the low-lying electronic branch would 
become strongly anisotropic (\B\ -C \A\): 

Hp) = - | ^ | -25c, + 3 -A(c x + c y ). (14') 

Returning to the full electronic Hamiltonian, let us recall that the Hund's term (|2|) is 
the largest one: we assume below Jh \A\,gQ. Therefore in all cases the electronic 
spectrum is shifted up or down by the energy ±J H S, depending on the e 2ff and t 2g mutual 
spin direction. At the ferromagnetic alignment of the all local and itinerant spins each of the 
two branches of Eq. ( |I"3"1) merely splits into two by adding the ± J# < S > energy term. For 
an antiferromagnetic (AF) ground state one needs to calculate the energy spectrum again. 
It is done below for the important case of two "canted" antiferromagnetic sublattices when a 
ferromagnetic spin component Mj coexists with the staggered magnetization, Sj = (— l) l Sj. 
We assume that antiferromagnetic ordering with the structure vector q, runs along the 
c-axis, q = ^(0, 0, 1). Such an AF - structure is often dubbed as the A-phase. 

We will perform our calculations treating local Mj and Sj as the classical vectors: (Mj + 
Sj) 2 = S 2 ; | S |= 3/2. The C2 g - electron wave function in the momentum representation 
acquires the form: 







1 a pa [ 




H 








[ Pp° J 



[ ^2 J [ /?p CT J 

where now a pa ,j3 p(T have a nontrivial spin-index dependence. Let S(q), the staggered mag- 
netization component, be directed along the z-direction, while M is parallel to the x - axis 
(in the spin space). In the calculations below we suppose q is along c-axis. The Brillouin 
zone is reduced along the c - direction, in accordance with a choice of q - vector. The energy 
spectrum is determined by the following equations: 
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Elp pa = (t(p) - J H a x M)lp pcT - J H ^zSlp p+ qa, 
Elpp+qa = (*(p + Q) - JH^xM)^ p+(l<J ~ JnCTzSlppa- 

Multiplying each equation by a z and eliminating ip p+Cia , one arrives to the following two- by- 
two matrix equation: 

E - t p+q - J h v x m] [E - t p + J H a x M] ^ p = J|5V P 

or 



E 2 - 2t Px , Py E + 2J H a x Mt Pz + P - P] ^ p = J 2 H S 2 %, 



(16) 



where we have separated in the t - matrix its dependence on the momenta parallel and 
perpendicular to the direction of the q - vector: 



and also we have used that i Pz+qz = — t Pz . Equation ([Tq) gives E = ±JhS in the zeroth 
approximation, as expected. Taking E = —JhS + e, we obtain: 



e 2 - 2J H Se + 2t Pvx>py (J H S - e) + 2J H a x Mt Pz + P pxjPy - P pz % = 0. 



(16') 



For the next approximation (e ~ t) in (|1 6 "| ) : 



iip p 



M. 



t PxiPy + cr x cos (9/2)i Pz ip 



p- 



(17) 



i.e. one obtains the generalized Anderson- Hasegawa term for the motion along the c - axis. 
Eq. ([17]) is the four-by-four matrix: two-by-two due to the orbital degeneracy and two-by- 
two due to the spin variables. After diagonalization over the spin variables, it provides the 
following tunneling matrix: 



£^ = ^±008(0/2)*,,. 



(18) 



To obtain the new spectrum one merely has to make a substitution in Eq. fll3|) : c z — > 
cos (|)c 2 (the number of energy branches is doubled because of reduction of the Brillouin zone 
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at doubling the periodicity along the c - direction. If, in addition, there is a homogeneous 
JT deformation, all t PxjPy and t Pz must be changed in the same manner as it has been 
done above). In accordance with the DE - mechanism, any connection between adjacent 
ferromagnetic layer in the A-phase fully disappears in the approximation (|T8|). To find 
further approximations, one has to take into account small terms of the order of t 2 / Jh <C t. 
The analysis below produces the scale for t in (|13|): | A |~ 0.16eV , i. e. these terms should be 
reasonably small and may be neglected except for some more peculiar phenomena in doped 
manganites. Regarding the magnitude of the JT-component, Eqs. ([|, |7|), it is suggested 
below that main physics to large extent depends on competition between the energy gained 
in the JT-distortions and kinetic (band) energy of electrons. Correspondingly, the JT-term 
is generally assumed to be of the same order as | A |, gQo ~| A |, in our model. 

IV. PARENT MANGANITE AS A BAND INSULATOR 

At low temperatures the stoichiometric ends of materials in series Ai_ x B x Mn03 are 
insulators. For the "right" end (x = 1) it causes no surprise: all manganese ions are Mn 4+ , 
i.e. comprise of the localized t 2g - spins. The magnetic properties are then described by an 
effective Heisenberg Hamiltonian commonly attributed to the superexchange mechanisms 
along the Mn-O-Mn bonds (the Hamiltonian itself may be constructed by the so-called 
"Anderson-Goodenough-Kanamori rules" (see in |^JT^JT^]). For the same approach to be 
applied for the "left" end (x = 0, LaMnOs), one should similary treat the Mn 3+ e2 g - shells as 
localized spin orbitals. All pertinent properties are often interpreted in terms of generalized 
microscopic Hubbard model (see e.g. [[n||15|]). The key feature of the Hubbard model is the 



assertion that for two electrons to be placed on the same site the energy cost is very high 
(the famous Hubbard U U" > due to the on-site Coulomb repulsion!). 

The Hubbard Hamiltonian approach has been challenged in ||. First, there are ex- 
perimental motivations for such challenge. For instance it was shown than in "doped" 
manganites, Lai^Sr^MnOs, at rather low concentrations, say, x = 0.2, the system displays 
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excellent low T metallic behavior Meanwhile, the nominal number of e2 9 -electrons per 



Mn site, N — 1 — x < 1, was changed only by one fifth in this study. From the theoretical 
side, it is of course clear that dealing with the strictly atomic d-orbitals would be a strong 
oversimplification. If a Mn ion is placed into the oxygen octahedron environment, the e2 g - 
terms are formed by the whole ligand, so that the "pure" (^-functions become considerably 
hybridized with the surrounding oxygen states (e.g. see the discussion in []10| , |T8| ). Hence 
the electronic polarization would undoubtly reduce the magnitude of the "Hubbard"- like 
(on-site) interactions. Finally, for the JT - effect (which itself is nothing but another form 
of the Coulomb interaction) to come up there is no need to use the Mn 3+ localized states 
picture as it will be seen from what is done below. The JT effect may persist even in the 
band model. 

Therefore, in what follows, we adopt the band approach to describe the ground state of 
LaMnOa. This approach rationalizes the major features of an insulating state in LaMnOs 
very well and merge into metallic state of "doped" manganites (see [^(J ) . 

Experimentally it is known [[I]] that LaMnOs transforms into the AF state below the Neel 
temperature, T N ~ 140K. Its magnetic structure corresponds to the A - phase, i.e., consists 
of ferromagnetically aligned layers which are coupled antiferromagntically along the one of 
main axes (that axis is chosen as c-axis below). As for the lattice, the cubic cell becomes 
compressed along the same direction because the oxygen octahedra exercise the collective 
JT ordering: the octahedra are elongated along x- or y-direction and are packed as it is 
shown schematically in Fig. 2. Correspondingly, the periodicity is doubled in the (x, y) 
plane by the JT - distortions, and is also doubled by the A - phase magnetic superstructure 
along the z-axis. Simple considerations will show that calculating the band spectrum for 
such a superstructure and filling the bands at one electron per manganese site immediately 
lead to the band insulator ground state. 

We preface these calculations with a brief discussion of the role of DE mechanism for a 
more simple model in which each "manganese" site has only one orbital for mobile electrons 
(in addition to localized spin Sj, of course). Instead of fll3f) , one has now the tight binding 
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spectrum: 

e(p) = t(cos (p x a) + cos (p y a) + cos (p z a)). (19) 

Imposing the Hund's term (^), Jh ^> t, the energy of an electron in the "half- metallic" (fully 
polarized) state, £& m (p), is: 

£(p) = -J H S + £(p) (19') 
and the electronic contribution to the total energy: 

E% = J eUp)n P j0^ = -JhS. (20) 

(per cubic unit cell). 

Rederiving Eq. fllTf ) in the simplified case described above for the A - phase, it is easy 
to find the electronic energy of the A - phase (M ~ in (TO)): 

^ = -I^rT4^^-^-Tts- (20,) 



From (|20|) and (|20'| ) one may conclude first, that terms of the order of t 2 /Jn in the 



electronic energy make the AF - state (the A - phase) energetically more favorable. That 
gain, being small, should be considered together with the superexchange contributions, if 
any. Moreover, both Eqs. (|20"D and (|20'|) were obtained for the "stoichiometric" case of 
one electron per Mn site. The DE mechanism in cooperation with large value of Jh makes 
all electrons have only one spin direction. Therefore, the bands in ( |20| ) and ( [2Q'| ) are fully 
occupied. In both cases (ferromagnegnetic and antiferromagnetic states) one would obtain 
the insulating phase for that model. 

Returning to the more realistic situation of the e2 9 -level, the two bands, and the A - type 
antiferromagnetic order, it follows again from Eq. (|18|) that in the main approximation the 
problem reduces to constructing the ground state of a single ferromagnetic layer. Unlike the 
previous case, a "simple" ferromagnetic state would be metallic because of the two bands: 



e±(p) = -\A\(c x + c y ) ± \A\yJ<? x + c 2 y - c x c y (13') 
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(since \B\ <C \A\ as was mentioned above, the -B-part in (|13|) was neglected) and one electron 
per cell. The Brillouin zone will be only half-filled (all electrons are polarized). Superstruc- 
ture imposed by the JT deformations, as shown in Fig. 1, suggests that since the Brillouin 
zone is now reduced by a factor of two, the same number of electrons may fill up the reduced 
zone, producing an insulator at the appropriate energy spectrum. Experimentally, LaMn03 
is an insulator and it is interesting to find the electronic band spectrum in the presence of 
the cooperative JT - distortions and then discuss the spectrum properties together with the 
experimental findings. 

Before we proceed further, let us make one remark. As it is well known, two dimensional 
ferromagnetic state is never stable being destroyed by spin fluctuations. Stabilization of the 
A - phase comes up due to the small terms in energy, which are responsible for remnant 
interactions between layers. Our estimates below for | A |~ 0.16eV and JjjS ~ 1.5eV would 
place the terms ~ t 2 / Jh on the scale of ~ lOOif which agrees well with the low value of the 
Neel temperature, T N ~ 150i^. 

However, "superexchange" terms, J • Sj, with J, being on the same scale as ~ t 2 / J H , 
may also become important. Role of these interaction at the small "electron" doping (i.e. 
close to the CaMn0 3 end) has recently been discussed in terms of the band picture in [[ZT|]. 

To calculate the electronic spectrum of our model in the presence of the "antiferrodis- 
torsive" JT collective deformations shown in Fig. 1, let us note that the structure vector 
of the superlattice is q = ^(1, 1, 0). The superstructure shows up in the JT - term (|7|) and 
depends on which of the two modes (Q2, Q3 or even their superposition) of Eq. (jfy is chosen 
to correspond to the JT local deformations on each of the two sublattices. Experimentally, 
the arrangement in Fig. 1 is close to the one in which octahedra are elongated along each 
x- or y- axis preserving the tetragonal symmetry in the perpendicular plane (For that 
mode the number of "short legs" Mn — O equals four and the "long" ones equals two. These 



deformations are also seen in some doped materials (e.g. see in |22|). Alternations of that 
specific type would be reflected in expression (]!]) as changes from 9 = 2n/3 to 9 = — 2%/3 
between the two sublattices. That will lead to the band's secular equation which can not 
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be solved in the analytic form. However, one can consider the contribution, which comes 
from the Q2 mode only; this deformation just changes the sign on the adjacent sites of two 
sublattices, as also follows from the Eq. (||) and Fig. 1. For that pattern one may easily 
right down Hjt as: 

(oA 



gQc 



exp (zqr) 



gQo 



exp (iqr)r x 



(21) 



V 10 / 



and obtain the spectrum by following the line of our reasoning used above, when we were 
deriving the Eqs. Ql6| , [17D (f x is a "pseudospin" Pauli matrix defined on the basis @). 
Writing explicitly 



e-t(p) V>p = r x V P+q , 



(by definition t(p + q) = — i(p)). Eliminating ^ p+q from Eqs. ([22]), one obtains: 



(22) 



£ T x + t(p) 



e f x - t(p) 



T x A 



The spectrum branches are then obtained from the zeroes of the two-by-two determinant: 



det 



e — 



e - tlpV + e 



t(p), f x 



(23) 



where e is the unit matrix, t(p) = t(p) f x , and 



*(p), f x 



i(p) ~ r x t(p)f x . 



(24) 



With i(p) expressed in terms of the "pseudospin" Pauli matrices: 

1 

*(p) = ~/+(p)e + ^f+(p)r x + —f-{p)T y , 



(25) 



where /+(p) =| A \ (c x + c y ), f-(p) —\ A \ (c x — c y ), after simple calculations, determinant 
(f23|) transforms to the form: 

, 2 - (sQo/2) 2 - f P + + \fl] + if f + f- P + - ieyfif. 



s 2 - {gQ Q /2f - \fl + f - iff + f_ 



0. 
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The resulting bi-quadratic equation produces the following four branches, £j(p)(z = 1..4). 
Each of these four branches is determined in the reduced Brillouin zone: 

biAp) = {(gQo/2) 2 + \f + + \f 2 ± ^f 2 [(gQ /2y + fl} + n} 1/2 

S3 A (P) = {(gQo/2) 2 + 5/4/ 2 + 3/4/ 2 ± ^3/ 2 [(sQo/2) 2 + ft] + fi)^ . (26) 

At large enough \gQo/2 | the branches £i,2(p) are not crossing two other branches, £3,4 (p). 
Filling them up by two polarized electrons per doubled unit cell, completes the proof that, 
indeed, insulating LaMnOs may be considered as band insulator. 

For example, two sets of the spectrum branches (Effi), £1,2 (p) and £3,4 (p) begin to overlap 
for £j=i(p) and £»=3(p) at p x = p y = tt/2. The overlap is direct which imposes some limit 
on the value of the JT mode which makes LaMnOs to be an insulator: 

I gQ \> 0.1 I A I . (27) 

The optical gap, hence, corresponds to excitation of an electron from the £j=i(2) (p) band 
into the £ i=3 (4)(p)- 

The optical spectra would, in principle, allow to access the gap directly. The detailed 
analysis of the optical absorption in the pure LaMn03 may be further complicated by the 
fact that formation of an electron- hole pair by light, strictly speaking, should not lead to 
creation of individual electron- and hole-like excitations due to inevitable polaronic effects 
in the ionic crystal LaMn03. Since the experimental level of current optical studies (see 
|PHl ) in LaMnOs can not address the issue yet, we shall not dwell upon further theoretical 



analysis. 

To complete this Section, one more simple comment may be helpful. Namely, while in 
the case of a single ion with one electron occupies the local degenerate e2 9 -level, it inevitably 
leads to the local instability caused by the linear JT term (|J), the ferromagnetic state with 
the band spectrum ([13]), filled up to some level, would remain stable with respect to small 
enough JT distortions. There is a threshold value for the magnitude of the JT - terms, 
before the distorted JT state may set in. This threshold is determined by the electronic 
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kinetic energy gain and the elastic lattice energy. In case of LaMnOs the existence of the 
cooperative JT deformations is confirmed experimentally. 



V. DOPED MANGANITES. PERCOLATION 

With the substitutional doping, Ax^B^MnOs, say in, Lai-^Sr^MnOs by Sr, the ma- 
terials initially show no metallic behavior with resistivity increasing by many orders of 
magnitude at T => 0, In metals, one would expect that resistivity behaves like 

p(T)=p + AT a , (28) 

where p , the residual resistivity, is due to the structure defects or impurities, while the 
T-dependence comes up from scattering on thermal phonons or from electron-electron inter- 
actions. In the second case, the electronic relaxation rate l/r ee ~ T 2 /Ep, while the phonon 
mechanism prevails at l/r p h ~ T 3 /9 2 > l/r ee , i.e. at T > 9(9 /Ep), where is the character- 
istic phonon frequency. The latter in manganites is about four hundred degrees Kelvin while, 
as we will see, for Ep in manganites one gets the scale Ep ~ O.leV. Hence, the electron 
mechanism might become important up to rather high temperatures. However, at elevated 
temperatures a lot of new effects, related to the colossal magnetoresistance phenomenon, 
start to play the dominating role. 

To understand the nature of the ground state of pure and doped manganites, we need, first 
of all, to concentrate our attention on the low temperature properties. It has been shown 
first in p8[ that low T resistivity changes its character with doping, so that the metallic 



behavior fl28|) sets in at x = x cr ~ 0.16 — 0.17 p^ , |25| . Interestingly enough, the onset of 
conducting regime coincides with the onset of the low temperature ferromagnetism, strongly 
indicating in favor of DE mechanism for the latter. 

The value of po may vary, even for nominally the same material and composition, im- 
plying that quality of samples may need further improvement. At the same time it is clear, 
that in all cases the change from insulating to metallic behavior takes place when x is close 
to the critical concentration, x cr . 

16 



The origin of that threshold and its value have been first understood in ||14}] in terms 



of percolation theory (for review on percolation theory see e.g. the book PB|). The latter 
considers any process which, roughly speaking, corresponds to some exchange between two 
adjacent local sites. A material, Ai.^B^MnOs, is commonly prepared by various methods 
at high temperatures. As the result, position of atom B, which is a substitution for a 
parent atom, is compelely random. Divalent atom B, locally creates a "hole" on adjacent 
Mn - sites which could serve as a charge carrier if there were no long range Coulomb forces 
in the dielectric phase, which keep the "hole" close to the negative charge at B~. When 
concentration is small, average distances between B atoms are large, which makes holes 
remain isolated. In the theory of percolation one may look for the concentration at which 
the nearest neighbor B - atoms start to form infinite clasters piercing the whole crystal. For 
the B - atoms on the cubic sites it is known as the "site" problem: on the simple cubic lattice 
it gives the critical concentration value at ~ 0.31. However, this is not exactly our case. For 
the doped holes, gathering on few manganese sites around an atom B, charge transfer takes 
place only along Mn-O-Mn bonds. Therefore the picture of a critical cluster, constructed 
from the B~-ions must be corrected: such cluster would already have a finite "thickness" 
due to the holes spread over surrounding Mn - sites. Numerical studies on the percolative 
models f27j have shown that this circumstance (i.e. the presence of the scale of a few lattice 
constants) strongly decreases the value for the critical concentration, rapidly converging to 
its value for the homogeneous problem, x cr ~ 0.16. The phenomenon of percolation in 
disordered media is akin to the other critical phenomena: it bears a singular behavior in the 
vicinity of x cr , \ x — x cr | -C x cr . Thus, it is expected that conductivity above x cr or Po 1 (x) 
in ([H]) should behave as a(x) oc(x — x cr ) 7 , where 7 is a critical index. There were not many 



attempts to verify singular behavior in conductivity. The data [p9].p8f are consistent with 



the index 7 ~ 0.5 - 0.6 [22 



The above arguments have been given at T = 0. One may try to apply similar ideas 
to the CMR phenomena at T = T c (see, e.g., in [§,0]). Indeed, if resistivity, p(T > T c ) is 
large enough, as it usually is, one may approximately take the conductivity: o~(T > T c ) = 
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p 1 (T > T c ) ~ 0. The fast increase in cr(T) at T < T c is then expected to correspond 
a(T < T c ) oc (T c — T) 7 . This is approximately true for La .8Sr .2MnO 3 near T c with 7 ~ 0.6 

The correctness of the percolation theory views may be also verified by independent 
measurements of a(T,x) and magnetization, M(T,x), or D sti ff(T,x) - the so-called "spin- 
stiffness" which determines the spectrum of the long-wave magnons: 

w(k) = D stiff k 2 (29) 

(D itself is proportional to M). The Kirkpatrick's relation gives: 

aocDM. (30) 

Eq. (|30D may be verified at x close to x cr or at T near T c . 

For percolation in the continuous phases description the physical picture is rather well 
understood. One of the helpful observations is that for the three-dimensional problem there is 
a concentration range in which percolation (infinite clusters) may take place simultaneously 
for two phases (insulating and conducting, for instance) though taking into account the 
surface tension at the boundary between two phases seems to impose limitations on relative 
thickness of the phases. 

In case of manganites, understanding on the microscopic level remains far from being 
complete, though the concept is by all means correct. In the above, the mixture of "two" 
phases, depending on the concentration, x, may look as intervened tiny "islands" and "lay- 
ers" of different "phases" , whose thickness even difficult to access quantitatively. Well below 
and well above the threshold concentration x cr , one may imagine each corresponding phase 
as a bulk formation into which the second phase is sparsely embedded. If there is a spill-over 
of charge carriers between two phases, it is the electro-neutrality condition which regulates 
the tiny domain sizes. 
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VI. PHYSICAL PROPERTIES IN TWO-BAND MODEL 



As it was explained in the previous Section, the percolation being a critical phenomenon, 
may cede soon to the onset of the homogeneous ferromagnetic phase, when screening becomes 
effective. In fact, some good samples of Lai-^Sr^MnOa show low temperature resistivity in 
the range 10~ 4 — 10~ 5 (Q — cm) ||45|| . Interest in applications of CMR was the main reason 
why the properties of manganites with x ~ 0.3 have been elaborated rather carefully 
(temperature of the resistivity peak, T p , reaches its maximum at 300 K for concentrations 
in this range). There are still indications [KJ that remnants of the second phase may still 
persist at x ~ 0.3. However, below we will analyze available low temperature experimental 
data for doped manganites from the band approach. Our conclusions will be that the band 
model (or Fermi liquid) approach seems to be valid at low temperatures, especially for the 
high quality samples. Data are hindered by sample's quality with local inhomogeneities or 
the second phase inclusions forming the scattering centers. For a number of compositions, 
depending on tolerance factor, conductivity is close to its value in the mobility edge regime. 
The two band model, thus, turns out to be a good starting point even at given doping level, 
providing together with the interpretation of the parent material, LaMn03, an avenue for 
unifying theoretical understanding of the manganites properties. 

Making use of the spectrum (|T3"D, we calculate the concentration dependence of the 
Fermi-level, Ep{x), the density of states (DOS), v(x), the spin stiffness, D(x), and the whole 
magnon spectrum, u(k,x), and the conductivity, a(u,x), both the Drude and the optical 
(interband) components. Theoretical results depend only on the single hopping integral, 
\A\. 

We start with the calculations of the spin wave spectrum. Let us write the deviations 
from the average spin, (S z ) for the localized ti g - spins (s = S — (S z )) as: 

S + (q) = (2(^)) 1 / 2 6(q), 5-(q) = (2(^)) 1 / 2 6 + (q), 

s z (q) = (S+S) q> (31) 
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(b + ,b - the magnon's operators). The first (8Ei) and second (SE2) order corrections to the 
ground state are calculated as perturbations in: 

V = -J H Si(at<TOi) = -J H SiHi. (32) 

i i 

For 8E2, the matrix elements in (R2|) are of the form: 



^p- k , T ;p| = -^(Tl s± ( k )U)x 

x («U< + P 1 ^; 1 ') , (33) 

(the coefficients a p ,/3 p are defined below). As for 8E\, its only role is to secure the proper 
behavior of the magnon spectrum at k — ► 0). One obtains: 

6E 2 = 2J 2 H Y / (S z )b + (k)b(k)x 

k 



I „M' _L Rl R*l' 1 2 



\a° a 



x 



Y.[Y. 'i7:r + J&r, i. (34) 



^ £{(p)-£[( P + k) 



where 



4|(P) = =F-/tf(&>+eM'(p). (35) 



Both sums in (|3^) run over Z, /' = ± (Eq.( [13])). The summation over I and p is limited by 
the occupied states (f ) only. The coefficients (a 1 , /3 p ) above are for the Bloch's states on 
the basis (0): 

< = (E^E^) 1 / 2 , ^' = ± (E^jE^) 1 / 2 , (36) 



(here E12, E21 are the off-diagonal elements of the hopping matrix £(p) in (|TT|) on this basis). 

Assuming in ( |3l| , |35| ) (recall that J H ^> |A| (JhS z is of the order of 1.5 eV |L| , while 
estimates below produce for \A\ ~ 0.1 eV), expanding (p 7 ^) in |A|/Jg would get a series of 
the Heisenberg spin Hamiltonians accounting for interactions with the increasing number 



of neighbors. (For a single band it was first noticed in |3(J; in this paper we are using the 



realistic two bands picture). After a somewhat tedious, but straightforward calculation, the 



first order term in | A | from (|34|) is ({S z ) = 3/2): 
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hw{k) = \A\(3 - c x -c y - c z )D{x)/3 



(37) 



and D(x) = D(E(x)) is given by the integral: 

d 3 p 



(27T) 5 



(here E is in units of |A|,pj = api). Quantum fluctuations may change the k - dependence 
in Eq. ([37]). 

Let us turn now to the calculation of conductivity, cx^u;, x), which is described by: 
e 2 h 2 



'j 



Voo 



£/o(k) [1 - /o(k')] ty(k) | v t | V(k'))^(k') | v 3 | ^(k)> 



x 



k,k' 



[5(e(k>) - c(k) - u,) - (e(kO " 6(k) + W )] 



(38) 



where /o(k) is the Fermi distribution function, k is a quasi momentum, %i is a velocity 
operator. For the cubic crystal = a^. 

We determine first the velocity operator, v = r (see ||31|| ): 

l&i(k) , i 



v(k) 



h dk h 



+ -[e,(k)-e,,(k)] (/k|ft|fk>. 



(39) 



The off-diagonal operator f2 is defined by the relation: 

(/k|^|/'k)=z| M f(r)^ 3 r 
and «k(r), the periodic Bloch functions on the basis of Eq. @ are: 



(40) 



exp|zk(aw — r)] 



x 



x I ot^j (r - na) + /^0 2 (r - na)j . 



(41) 



With the one-site integrals only in ( flOD and Eqs. (|36|) : 

<Zk|«|rk> = (-sirx ^)( C3/ - 



h 4 



M 2 



(42) 



Matrix elements in (|39| , [42]) produce transitions from occupied parts of the £+(p)-band into 
empty states in the e_ (p)-band. 
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With all the above, we arrive to the Drude (intraband) contribution which in the clean 
limit is: 

e 2 \A\ 

o-Drude(u,x) = 2ir——j5(uj)I Dr (x), (43) 
San 

IdAx) ~E / <|V p£ (p)|, (44) 

(the integral in lDr{x) is over the Fermi surfaces). 
The "optical" (interband) contribution is 

3vre 2 1 r d 3 p . 

x n(e+(p))[ 1 - n(e_(p))] • 5 {Co - 2R(p)) , (45) 



a opt[ u, x) = — ^ y j-y 3 sm M c, - c,J x 



where \A\Cb = %u (Eqs. ( f4"4| ) and ( ^5| ) agree with the similar calculations of |S2|). For the 
low temperature spectral weight 

oo 

iV e// = —a 3 / a(cj)rfw (46) 
o 

one obtains both the Drude and the interband contributions, respectively: 



= ^rl^l^W. JV.T/ = ^\A\%t(x), (47) 



(with Iopt(x) directly obtained from fl45|) ). In Fig. 2 we plotted our results for the Fermi level, 
Ep(x) = \A\E(x),v(x) = u(x)\A\, D(x) and Idt{x) (the shaded area shows schematically 
the concentration range for a percolative regime). In Fig. 3 evolution of the Fermi surface 
with x is shown (Fermi surfaces in Fig. 3a, 3b, 3c are for concentrations x = 0.2, 0.3, 0.5 
respectively. An interesting fact is that x = 0.3 is the concentration at which a "neck" 
develops at the zone boundary of the Fermi surface. In other words, concentration x = 0.3 
is the point of "2.5" - Lifshitz transition at which the change in the Fermi surface topology 



occurs. 



VII. THE MODEL AND EXPERIMENT. DISCUSSION 

In the discussion below we restrict ourselves mainly by the "isotropic" materials, i.e. 
by the range of concentrations x < 0.4. Around x ~ 0.5 new phenomena start to de- 
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velop in doped manganites, such as charge ordering fl33U34fl, A- phase in Pr .5Sr .5MnO3 and 



Ndo.45Sr .5 5 Mn03 p5|j36[1 , structural transitions caused by magnetic field in Pro.5Sro.5MnO 



3) 



Ndo.5Sr .5Mn0 3 |37],|38],(4T| and "spin valve" effect in Nd .46Sr .54MnO 3 Though some 



preparatory steps have been made to address these issues, e.g. Eqs. ( [H| , [H]) in Sec. Ill, 
the purpose of this paper is mainly to access the limits of applicability of the band structure 
theory above, and its parameters. For these reasons the concentration range x ~ 0.3 is 
chosen where more data are available. A word of caution still would be helpful. Although 
the quality of materials has been dramatically improved for the past couple of years and 
many results are becoming quite reproducible, systematic studies of the same phenomenon 
on the same materials, which would have been performed in different laboratories, are rather 
rare. It is also not well known how data, obtained on single crystals may differ from the 
ones obtained in crystalline films. 

The spectrum (|3|) and the results above, such as (|37j , [43| - f47D include the single en- 



ergy scale, | A |. We choose to determine this parameter from the measured spin stiffness 
coefficient in (|3"TD: 

Mk) = I A I D(x), (38') 

6 

1. e. 

a 2 

Dsuff = j I A | D(x). 

Being a long-wave characteristics, a latter must be less sensitive to defects of sample's quality. 
We use a ~ 3.86A, D(x ~ 0.3) ~ 0.45 (according to our results in Fig. 2) and the data, 
collected in |I2] (Table 1). The results for | A | show that the bandwidths for different 
materials, W = 6 | A |, do not vary significantly (W ~ 0.7 — l.OeV). This is in favor of the 
fact that the tolerance factor itself which is responsible for variations in the Mn — O — Mn 
bonds angle for the different materials, is not of much importance, as it has already been 
suggested for the bandwidth in P0] , P2]| . As for the concentration dependence, D(x) in Fig. 

2, it reflects the fact that number of electrons and, hence, the saturation moment M(x) 
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decreases with the increase of x. There are measurements of D s uff(x) [0,0 which show 
some pronounced dependence on x < 0.2 and saturation at x — 0.28,0.3 |43|]. The decrease 
in D(x), according to our results, shown in Fig. 2, at x ~ 0.3 is below the experimental 



accuracy in [4"3]]. As for measured dependence D st iff(x) for smaller concentrations, here one 
may meet the range for the percolative behavior, shown in our Fig. 2 by the shaded area. 
We will return to it later, although it is worth mentioning that data |43|]44]] do not agree 
with each other above x = 0.2. 

With | A |~ 0.16eV (Lai-xSr^MnOs) one may find the thermodynamical parameters 
and the density of states, u(x) = v{x)j \ A |. Correspondingly, the Sommerfeld coefficient, 
7, for our spectrum (|i~3"l) is 

7 = tt 2 z>(x)/3 \A\. (48) 

An attempt to investigate properties of Lai_ x Sr x Mn03 systematically has been done by 
Okuda, where special efforts were applied to single out the electronic component into the 
specific heat by subtracting both phonon and magnon contributions. According to |23| , 7 ~ 



3.5mJ/mole ■ K 2 , while v(x) ~ 0.45 at x = 0.3 from Fig. 2 produces 7 ~ 6.2mJ/mole ■ K 2 . 
Note, however, that v(x) has a kink at x = 0.3. That kink takes its origin from the fact 
that this concentration is the point of the Lifshitz singularity, as it may be clearly seen from 
the Fermi surface pattern in Fig. 3. The total magnon contribution into specific heat is 
proportional to T 3//2 . The proximity to the Lifshitz "2.5" transition results in appearance 
of the same T— dependence of the electronic specific heat. As the result, the procedure, 
24" I, of extracting magnon T 3//2 terms becomes less transparent. It is worth mentioning 



that unambiguous determination of the electronic contribution is known |T[ as a difficult 
task. The low T electronic specific heat of a few other manganite compounds has also 
been measured with 7 in the 3 — 8mJ/mole ■ K 2 range. Taking into account the possible 
complications with the rapid energy dependence in the density of states close to x — 0.3, 
the accuracy of Eq. ( f4"8"D seems to us rather satisfactory. 

Important issue is a sensitivity of all data to disorder, which is inevitably present at 
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the substitutional doping. The concentrations are not small by all means. The residual 
resistances in Eq. ( p8|) obtained for nominally the same compositions, may vary significantly 
for different groups of compounds. So, it is clear, that to some extent, the best values of po(x) 
are still to be determined. Meanwhile, it remains unclear whether po{x) only characterizes 
the sample's quality or there is an intrinsic component in the residual resistance. In favor 
of the former suggestion tells the results of recent experiments performed by Hl9| , ^5| , ^6| . 
New records have recently been reported for conductivity, measured in a few compounds. 
Among them there are resistivity data for Lai_ x Sr x Mn03 in crystalline films which give the 
value for po(x) as low as 10~ 5 f2 • cm |[45|| . This is a typical metallic conductivity range. 
Making a substitution in Eq. (|43|) for Ti5{u))— > r/(l + (cur) 2 ), we obtain in this case 
h/r | A |~ 2-10 -2 , i.e. for the inverse life time, h/r ~ 3-10 _3 eV^(~ 30i^). At the same time 



the residual resistivity for other materials studied in [45], lie in the range 100 -r-300(/ifi — cm), 
i.e. h/r | A |~ 0.5. These findings become more transparent being expressed in terms of 
the mean free path. With the average velocity of an electron on the Fermi surface 

v =< v 2 > 1/2 = (| A | a/K)(2I Dr (x)/D(x)) 112 , (49) 

the mean free path, I = vr, is typically ~ 3a for materials with p ~ 10 _4 (fi • cm), while in 
the best LSMO samples it is around 80a. Whether the values of po mentioned above may be 
or may be not improved by a more careful sample preparation process, remains to be seen. 
In any case, for some current materials |45|] the conductivity regime lies close to the mobility 



edge. If the short mean free path in these materials is an intrinsic feature, it may be related to 
the local fluctuations in the tolerance factor (JJ) caused by the difference in ionic radii at the 
«SV— substitution. Indeed, with < r >= 0.14nm, < r$ r >= 0.129nm, < r^ a >= 0.136nm 
for Pr 3+ and Nd 3+ one has < r Pr >= 0.129nm and < r Nd >= 0.124nm Jl[(p. 194). As it 
has been emphasized already, our two-band approach is a simple way of the Fermi liquid 
description, when interactions are not supposed to be remarkably strong. The strength of 
interactions is provided by the T 2 — term in resistivity, in which electron interactions come 
from Ti/r ee tr = X' (h/T ee ), where r ee is the total quasiparticle relaxation time and A' < 1 gives 
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the fraction of Umklapp processes. We use [ 47 1 : 



h/r ee ~ Xw 3 u(x)T 2 . (50) 



In ([5(]) A is a value of the interaction strength in terms of Ep. Using f4T^j4"6|| , one obtains 



for LSMO: AA' ~ 0.3, typical of good metals. As for the two other materials in [45|, the 
T-variations of resistivity scale in the magnitude with their residual resistivity and most 
probably are caused by defects 

T- dependence in the optical conductivity, cr(u), attracted recently much attention 
2^ , |45|^9| as a manifestation of changes in the conductivity mechanism from metallic to 



localized polarons at elevated temperatures. We discuss only a few results pertinent to the 
low T band mechanisms. First note, that the temperature dependence in u{uj) at T < 100 
K for uj around leV is most pronounced below leV (see Fig. 2 in j|5] ). This agrees well 
with our estimates for the bandwidths, W < leV. As for a quantitative analysis, there are 



problems of an experimental character. Assuming N e jf(u) [15] would give our N e f^s in Eq. 
([IT]) at uj ~ leV and that both the Drude and optical contributions are approximately equal 
we obtain N e ff ~ 0.25 which is reasonably close to the values in jKJ , lesser than N e ff 



for the single crystal data, La .67Cao.33Mn0 3 |4T| , and a factor ten bigger than N e ff for 



Lao.ySro.sMnOs in [[nj . We believe that such a difference originates from poor data for the 
optical <j(uj) in the "Drude-tail" range. 



Finally, the optical gap for pure LaMn03 was identified [|Tl|] at A ~ 1.2eV. A rough 
estimate from the band insulator picture (|26|) gives gQo/2 ~ 0.6eV, i.e. the JT— coupling 
is rather strong (| A \ ~0.l6eV). 

We conclude the discussion by a few comments regarding the spin wave spectrum. In 
La . 7 Pbo.3Mn0 3 [0 the spectrum fits well (|7J). Eq.(gTJ) follows from (|3§ at A < J H (S Z ), 
with quantum corrections neglected (we have also neglected terms of the order of t 2 / Jh)- 
Meanwhile, strong deviations from ( |3TD have been observed at £ > 0.25 along the (0,0,£)- 
direction in Pr .63Sr .37MnO 3 The spin stiffness changes only slightly from material 

to material [50|, and other low temperature characteristics including strength of electron 
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interactions, also seem not to vary much. Unlike [5(J, we suggest with |H[ that differences 



in low-T spin dynamics for two materials do not correlate with their behavior at higher 
temperatures, but with a tendency to a low temperature charge ordering or formation of the 
A - phase at x close to 0.5 (we will address this issue elsewhere). 

Let us now return to the problem of how doped manganites proceed from the insulating 
end, LaMnOs, with the increase of concentration, x. At x < x cr doped holes remain localized 
in the vicinity of corresponding clusters. If the clusters are small enough they might be called 
"ferromagnetic polarons" but this is a semantics. In [deGennes, 1960] it has been argued that 
such isolated microscopic ferromagnetic inclusions do produce a macroscopic consequences, 
forming "canted" ferromagnetic moment even in the A— phase of parent LaMnOs. In other 
words, the effective radius of ferromagnetic interactions is larger than short range coupling 



in the percolative scenario for conductivity. The relation (|30|) based on the assumption 
of equal effective radii both for conductivity hopping and ferromagnetic interactions, is, 
thus, not expected to be applicable near x cr = 0.16: magnetic ordering and ferromagnetic 



component in the spin waves, may develope earlier. The experimental results i^jli ], indeed, 
show the stiffness coefficient, D sti ff(x) nonzero below x cr and starting to increase almost 
at x < 0.15. However, the further sharp increase above x cr = 0.16 with a saturation near 
x ~ 0.3, is the bright demonstration in favor of the notion of percolation: the increase in 
D s tiff{x) just follows increase in the amount of the ferromagnetic component. Note, that 
D S uff (x) is proportional to the value of magnetization, M(x). In the band picture, M(x) is 
(4 — x)hb- As it is known since the old results [||] that M(x) actually increases at x > 0.1. 
Interesting thing is that M(x) has dependence, close to (4 — x)fig, even at x ~ 0.3 — 0.4, 
indicating that admixture of "insulating phase" still persists at these concentrations. 

Another proof comes from recent pulsed neutron experiments f22||5"2|j . These experiments 



probe the local arrangement of the oxygen octahedra. With the random disorder, one may 
expect that the degeneracy of the e2 9 — terms would also be lifted randomly, providing a 
distribution in the JT deformations of the oxygen octahedra. Meanwhile, it has been shown 



in p2| that the pair distribution function displays the well pronounced peaks at the values 

27 



of Mn — bond which are characteristic for the parameters of the elongated octahedra in 
the parent LaMn03. According to p2"| , in La1_2.Srj.MnO3 the presence of such insulating 
inclusions is seen up to x = 0.35. The local neutron probe P_| , |__| can not resolve the size and 



the pattern of the distorted areas. One possible view may be to consider them as randomly 
positioned centers or the JT "polarons" with a scale ~ 15 — 20 A. However, we have to 
keep in mind, that in the 3D percolative regime, the percolative paths (infinitely connected 
clusters) may coexist for both phases simultaneously. The analysis of low temperature 
metallic properties performed in the previous Section, seems to indicate that at least for 
concentrations x ~ 0.3, the ferromagnetic phase component, in the first approximation, 
occupies most of the bulk with other phase embedded into it and also seen as reasonably 
small scattering centers. 

The Kirkpatrick relation (|30|) is expected to be approximately correct in homogeneous 
media regime, near T c , where onset of ferromagnetism may be interpreted as the simultaneous 
onset of a new (metallic) mechanism for conductivity. As it was pointed out at the end of 
Section VI, Eq. ( |3"0"D may be satisfied reasonably well. At the same time, it has been already 
known that in La .8Sr 2 Mn0 3 the peak temperature, T p , and the Curie temperature, T c , 
do not coincide (T p = and T c =). These experiments were mostly considered in the literature 
as a crossover || in which conductivity regime via the thermal hopping of polarons localized 
by the thermal lattice disorder cedes sharply to a metallic regime with a short mean free path. 
In recent STM experiments [|53| it has been successfully demonstrated that even around T c 



the regime in Lao.7Cao.3Mn03 is better described as a coexistence of the two percolative 
phases: metallic and insulating ones. The characteristic scale for sizes of subphases having 
the shape of smooth clouds are from tens to hundreds of nanometers. 

Another experimental support comes from the Mossbauer spectroscopy measurements 



54fl . The strong paramagnetic signal has been observed at T > T c . Decrease in T leads to 
the appearance at T < T c of the strong signature of the ferromagnetic state, namely the "6 
picks" structure in the Mossbauer signal. However, even at T < T c the paramagnetic signal 
persists up to T ~ 20 K, which is much lower then T c . These studies strongly support the 
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exsistance of the two phases and therefore the percolative scenario. 



VIII. CONCLUSIONS 

We have discussed above two major issues in doped manganites: 1. - the concept of 
percolation, leading, as it was first predicted in ||, to the new understanding of the phe- 
nomena related to the phase separation, which takes place in these materials at low enough 
concentrations^. - the possibility to interpret some low temperature data in the habitual 
band (or Fermi liquid) picture. The latter makes possible to suggest an unifying view which 
circumvent the often raised concerns regarding the role of strong interactions, importance of 



which is commonly expected in the transition metal oxides [17]]. The Fermi liquid approach 
is a selfconsistent theory, which of course, may be broken at strong interactions. The only 
way to justify the Fermi liquid approach is to see whether the experimental results do fit the 
theory predictions or not. As it was shown above, the theory works reasonably well at x = 
and at x ~ 0.3 for Lai_ x Sr x MnO^, probably because the DE— mechanism results in rather 
effective de-localization while the JT - effect partially takes care of the Coulomb interaction. 
In the intermediate regime with respect to x, the tendency towards the phase segregation 
turns out to have a percolative character and, therefore, theoretical analysis becomes much 
more complicated. 

The question whether the samples quality is an intrinsic property or may be further im- 
proved, seems to be of principle importance. Experience with LSMO materials |4l| implies 



that attempts made to improve the crystals quality are worth of trying. Already the LSMO 
samples with 1/r ~ 20K are encouraging, being reasonably close to the possibility of study- 
ing metallic manganites by methods traditional to the normal physics of metals, including 
the possibility of the de Haas - van Alphen studies in fields of the order of 30 — 40 Tesla. 
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Figure Captions 

Fig. 1 In-plane staggered distortions inside the ferromagnetic layer. Solid segmenets rep- 
resent the elongations of octahedra. The new unit cell is shown as the dashed square. 
Fig. 2 The Fermi level, E(x), DOS, u(x), the spin stiffness coefficient D(x) and the Drude 
conductivity, l£> r {x), plotted as a function of concentration, x, for the spectrum, given in 
|3 . The shaded area shows the percolative regime where Eqs. (37-46) are not applicable. 



Fig. 3 The Fermi surface a) x=0.2; b) x=0.3; c) x=0.5. 



35 




0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 



X 



